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Abstract 

We study the resonances of the quantum kicked rotor subjected to an excitation that 
fohows a deterministic time-dependent prescription. For the primary resonances we 
find an analytical relation between the long-time behavior of the standard devia- 
tion and the external kick strength. For the secondary resonances we obtain essen- 
tially the same result numerically. Selecting the time sequence of the kick allows to 
obtain a variety of asymptotic wave-function spreadings: super-ballistic, ballistic, 
sub-ballistic, diffusive, sub-diffusive and localized. 
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1 Introduction 



The kicked rotor model is a paradigm of periodically driven systems with 
classical chaotic behavior and its quantum version, the quantum kicked rotor 
(QKR), can be considered a cornerstone in the study of chaos at the quantum 
level [I]. This topic is a matter of permanent attention, both theoretical and 
experimental, for example it has been recently reported [2] the experimen- 
tal realization of the quantum kicked top (a variant of the QKR.) to study 
signatures of quantum chaos. 



The behavior of the QKR has two characteristic modalities: dynamical local- 
ization and ballistic spreading of the variance in resonance [3]. These behaviors 
are quite different and have no classical analog. They depend on the relation- 
ship between the characteristic time of the free rotor and the period associated 
to the kick. When the dimensionless period of the kick, T, is an irrational mul- 
tiple of 27T the average energy of the system grows in a diffusive manner for 
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a short time and afterwards dynamical localization appears. In this case, the 
angular momentum distribution has a maximum at its mean value and de- 
cays exponentially. For T rational multiple of 27r the behavior of the system 
is resonant with ballistic spreading, and the angular momentum distribution 
evolves from the initial one in such away that its standard deviation has the 
time dependence a{t) ~ t. The two types of values of T determine the spectral 
properties of the Hamiltonian, for irrational multiples the energy spectrum is 
purely discrete and for rational multiples it contains a continuous part. Both 
quantum resonance and dynamical localization can be seen as interference 
phenomena, the first being a constructive interference effect and the second a 
destructive one. 

Different variants of the QKR have been studied, but with the aim of this 
work in mind, it is interesting to highlight the one developed by P|f5] where 
the QKR is subjected to two unitary operators with two different values of 
the strength parameter K. The dynamics of the QKR kicked according to a 
Fibonacci sequence is studied in Ref. This is done outside the resonant 
regime, where a sub-diffusive behavior is found for small kicking strengths. 
In Ref. [S] we have also studied the QKR excited with the same Fibonacci 
prescription but in the resonant regime. It was shown that the standard de- 
viation of the primary resonances maintains the well-known ballistic behavior 
but the secondary resonances have sub-ballistic behavior, characterized by the 
time dependence of its standard deviation a{t) ~ P, with 1/2 < 7 < 1. This 
change of behavior is similar to that obtained for the quantum walk (QW) 
subjected to an aperiodic Fibonacci excitation [B] where the ballistic behavior 
is lost and the sub-ballistic behavior appears. We have also investigated the 
QKR in resonance subjected to decoherence: a) with a Levy waiting-time dis- 
tribution [7|l8] and b) with an unitary noise defined through Kraus operators 
[H]. In both cases we find that the secondary resonances have a sub-ballistic 
behavior, while the principal resonances are not affected. On the other hand, in 
Refs. pUfTT] the QKR subjected to noises with a Levy distribution out of res- 
onance was investigated, showing that the decoherence leads to a sub-diffusive 
regime (7 < 1/2) for a short time before localization appears. 

The quantum resonances and the dynamical localization have been experi- 
mentally observed in Refs. [T^13fl4|[T5] . The QKR was experimentally real- 
ized through a dilute sample of ultra-cold atoms exposed to a one-dimensional 
spatially periodic optical potential that is pulsed on periodically in time (to 
approximate a series of delta function kicks). However, the ballistic growth of 
the secondary resonances was not confirmed. It was found that the momen- 
tum distribution profiles saturated to a final distribution but different from 
the characteristic exponential profiles of the dynamical localization. In Ref. 
[13] the usual QKR is modified, by adding other elements such as a quasi- 
momentum, in order to explain the absence of the ballistic growth. However 
later experiments [13] has confirmed this growth. 
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In a recent paper [T7j we have studied the QW model using a time-dependent 
unitary coin operator. We found different types of anomalous diffusion de- 
pending on the time prescription for the coin operator. In the same line of 
thought, in this paper we develop a simple QKR model with a generalized 
time strength parameter K{t) that allows an analytical treatment of the pri- 
mary resonances. We study the connection between the time evolution of K{t) 
and the type of spreading of the system. We analytically show how to select 
the time dependence of K{t) to obtain a predetermined power-law distribu- 
tion of the standard deviation. We propose this toy model to reproduce the 
experimental behavior reported in Refs. [TSfT^ . 

The paper is organized as follows: In the next section we present a brief re- 
vision of the QKR equations, in the third section we obtain analytically the 
standard deviation for the primary resonances, in section IV we choose a 
time dependence for K{t), in section V numerical results are presented for the 
secondary resonances, and in the last section we draw the conclusions. 



2 Quantum kicked rotor 



The QKR is one of the simplest and best investigated quantum models whose 
classical counterpart displays chaos. In this paper we use a modified version 
of the usual QKR where a deterministic time dependence for the strength 
parameter K[t) is introduced. The Hamiltonian has then the following shape 

p2 oo 

H{t) = ^ + m cos eY^6{t- nT), (1) 

n=l 

where the external kicks occur at times t = nT with integer n, T being the 
kick period (the index n will be equivalent to time in units of T, given the 
delta function appearing in Eq.(IT])), / is the moment of inertia of the rotor, 
P the angular momentum operator and 6 the angular position operator. The 
last two operators satisfy the canonical commutation rule 

[P, 6] = -ih. (2) 

In the angular momentum representation, P\l) = lh\l), the wave- vector is 

oo 

mn))= E «K^)IO (3) 

(=— oo 

and the average energy is 

oo 

Ein) = {^\H\^)=e ^ f \aiin)\\ (4) 

Z=— oo 
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where e = h'^/21. Using the Schrodinger equation the quantum map is readily 
obtained from the Hamiltonian Eq.(IT]) 

oo 

ai{n + l)= J2 Uijaj{n), (5) 

j=-oo 

where the matrix element of the time step evolution operator U{k) is 

Ui,=i-^^-'^e-'^'^^l'j^.i{K{n)), (6) 

Ji being the /th order cylindrical Bessel function and its argument the dimen- 
sionless kick strength n{n) = K{n)/h. With the aim to generate the dynamics 
of the system we will consider different time dependencies of the strength pa- 
rameter and combine the corresponding time-step operators U {K{n)) in a large 
sequence. The resonance condition does not depend on k and takes place when 
the frequency of the driving force is commensurable with the frequencies of 
the free rotor. Inspection of Eq. ([6]) shows that the resonant values of the scale 
parameter r = eT/2h are the set of the rational multiples of 27r, r = 27r p/q. 
In what follows we assume, that the resonance condition is satisfied, therefore 
the evolution operator depends on k, p and q. We call a resonance primary 
when p/q is an integer and secondary when it is not. 



3 Primary resonances 



We shall first consider the primary resonances. Using the recursion relation 
satisfied by the Bessel functions Eq.(|5]) is solved. Its general solution is written 
as 

oo 

Mn)= E HY^'^M Ji-An*), (7) 

j=-oo 

where aj(0) are the initial amplitudes and 

n 

n* =J2 '^M- (8) 

m=l 

From Eq.Q it is clear that the propagation speed of the probability ampli- 
tudes is given by An* /An = K{n). The probability distribution at time n is 
given by 

Pi{n) = \ai{n)\^, (9) 

that can be expressed as 

oo 

Piin)= E H)''''akiO)a*iO)Ji_,in*)Ji_jin*). (10) 

j,k=—oo 
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Using the properties of the Bessel functions all the moments of Pi{n) can 
be calculated analytically from Eq.f llOp . If symmetrical initial conditions are 
taken, the odd moments vanish due to the symmetrical dependence on the 
initial conditions. We calculate the first and second moments, obtaining: 



Mi(n) = -n* ^ 53 \aj{0)a*_^{0)] + Mi(0), 



M2(n) 



[n 



*\2 



1- Y: 3^[%(0)a*+2(0) 



j = -oo 



+n* (2j + l)53[a,(0)a*+i(0) 

j=-oo 

+M2(0), 



(12) 



where 3? [x] and [x] are the real and imaginary part of x respectively. The 



standard deviation cr = y M2 — Mf has the following long-time limit 

n 

(T{n) ^Cn* = CJ2 (13) 

m=l 

where C is a constant that only depends on the initial conditions. 



4 Choosing K{t) 



Among the alternatives to choose the time dependence of K{n) we study the 



case: 



(14) 



where a is a real parameter and is a constant. Substituting Eq. (fT4|) in 
Eq. (fT3l) and omitting a multiplicative constant the following asymptotic re- 
sults are obtained 



Inn, if a = 1, 



am — )■ < 



(15) 



n"^, if a 7^ 1, 
where 

7 = 1 — a. (16) 

The asymptotic spreading of the wave function shows seven different behaviors 
that can be characterized by the exponent 7: 
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(a) 7 > 1, super-ballistic. 

(b) 7 = 1, ballistic. 

(c) 1 > 7 > 1/2, sub-ballistic. 

(d) 7 = 1/2, diffusive. 

(e) 1/2 > 7 > 0, sub-diffusive. 

(f) 7 = 0, logarithmic. 

(g) 7 < 0, Bessel-localized. 

The result for 7 = 1 is expected because the system reduces to the usual QKR 
in resonance. In the other cases the system shows unexpected behaviors. For 
7 < the system does not spread when t — )■ 00, it becomes localized in 
angular momentum. We have called this behavior 'Bessel-localized' because 
the distribution decay is slower than the case of exponential localization, as can 
be seen from Eq. (fTO|) . Note that this equation with initial conditions afc(O) = 
Skm is reduced to: 

Pi{n) = [Ji_,^{n*)]\ (17) 

In the other extreme for 7 > 1 the spreading of the wave-function is super- 
ballistic. In this case not only the external frequency matches the natural 
frequency of the free rotator (resonance condition) but the intensity of the 
external strength grows with time as well. We could think that this condition 
produces a new kind of reinforced resonance. In order to visualize the asymp- 
totic behavior of the distribution Pi, we use Eq. (IT7|) to calculate the profiles 
of the distribution. These profiles are shown in Fig. [1] for several values of a 
after a time long enough has elapsed, to see the spreading of the distribution. 
We take the initial conditions ai{0) = Sio, to produce a symmetrical evolution. 
From this figure, it is clear that as a grows the spreading speed decreases, the 
distribution shrinks and the two extreme peaks come closer. When a = 1/2 
the standard deviation spreads out as a{t) ~ t^/^, but the probability distri- 
bution Pi is not Gaussian. This confirms that the evolution corresponds to 
a coherent unitary process and not to the classical diffusive behavior. In the 
extreme case when a > 1 the distribution is restricted to a small interval 
around its initial value and the two extreme peaks melt into a very narrow 
peak. It is important to point out that the distribution profiles presented here 
are qualitatively similar to those obtained in Ref. [TGlTTj for the QW. 

To conclude this section, it is important to point out that the specific form 
of the kick sequences in Eq. (fT^ is not essential to obtain the qualitative 
behaviors described above, they are independent of this choice. Note that cr(n) 
is determined by Eq. (fT3|) and the qualitative behavior of K,{n). If K,{n) — )■ 0, 
(j{n) cannot grow faster than in the ballistic case, therefore it behaves as in 
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Fig. 1. The probability distribution Pi, for a primary resonance as a function of the 
dimensionless momentum / at time n = 300 with the initial condition ai{0) = 6io. 
The values of the parameter a are, from top to bottom, —0.1, 0, 0.5, 1 and 2. 

the cases denoted above by (b), (c), (d), (e) or (g), depending on the rate of 
decay of /t(n). When K,{n) but remains bounded, the behavior is balHstic. 
Finally when n{n) is not bounded, grows super-ballistically as in case 

(a). 



5 Secondary resonances 

Now we inquire about the secondary resonances of the QKR. In the gen- 
eral case the analytical development is very cumbersome because the matrix 
elements of the unitary map in Eq.(l6]) include the time dependent phases 
g-«i eT/h^ These phases impede the integrability of the problem. It could be 
possible to obtain analytical results for some special function K{t) following 
the work of Ref . \)3\ ■ However, we opt to study the secondary resonances nu- 
merically iterating Eq.([S]) with the time dependent kick proposed in Eq. fll4l) 
to calculate the distribution Pi and the standard deviation. In Fig. [2] the stan- 
dard deviation a for the secondary resonance p/q = 2/5 is presented for several 
values of the parameter a. It is seen that a has a power-law growth with an 
exponent 7 that depends on a, but now the relation given by Eq. (fT6!l does 
not hold. The values of 7 were adjusted for the last hundred values of n. The 
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Fig. 2. The standard deviation (j{n) as a function of the dimensionless time n in 
log- log scales for the secondary resonance p/q = 2/5. The initial condition for the 
amplitudes is 0^(0) = Siq. For large n the curves satisfy, on the average, a power law 
rf. The values of 7 are, from top to bottom, (1) 7 ~ 1.4 for a = —0.1, (2) 7 = 1 
for Q = 0, (3) 7 ~ 0.5 for a = 0.1, and (4) 7 ~ for a = 0.2. 

ranges of a for the different behaviors are now the following: 

(a) a < 0, super-ballistic (7 > 1), 

(b) a = 0, ballistic (7 = 1), 

(c) < a < 0.1, sub-ballistic (1 > 7 > 0.5), 

(d) a ~ 0.1, diffusive (7 ~ 0.5), 

(e) 0.1 < a < 0.2, sub-diffusive 0.5 > 7 > 0, 

(f) a > 0.2, localized (7 < 0). 

If we compare these ranges with the corresponding ranges of a for the primary 
resonances we find a similar qualitative behavior, as in the cases (a), (b), (c), 
(d) and (e). However in the present calculation, the logarithmic case is not 
present and it is not easy to classify the type of localization obtained. 

In Fig. [3]we present the distribution Pi as a function of the angular momentum 
for four values of a. The distributions show a symmetrical spreading around 
a central peak. The spreading shrinks as a grows, going towards the localized 
behavior, as in the case of the primary resonances. 
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Fig. 3. The probability distribution Pi, for the secondary resonance p/q = 2/5, as a 
function of the dimensionless momentum I at time n = 300 for the initial condition 
a/(0) = 6iQ. The vertical scale is logarithmic. 

6 Discussion and conclusion 

The resonances of the QKR, as well as the dynamical localization, were first 
established numerically |lj and the analytical treatment of the resonances was 
developed by Izrailev and Shepelyansky [18]. Thirteen years ago the primary 
resonances were experimentally observed by Bharucha et al. pUl], in samples 
of cold atoms interacting with a far-detuned standing wave of laser light. 
More recently Kanem et al. [15] have observed secondary resonances. These 
notable series of experiments, together with those in Refs. [T2|14fl9j . have 
drawn much attention because they have allowed to study fundamental aspects 
of the interaction between matter and radiation, and may also contribute 
to the experimental basis of quantum computation. In Ref. [13] the usual 
QKR model was used to describe the experimental results where dynamical 
localization was found out of resonance. However, they did not observe the 
ballistic growth in resonance predicted by the model. Instead, an unexpected 
decay of the distribution was obtained. The variety of anomalous behaviors 
of the resonances presented in this paper could be useful to interpret these 
experimental results. 

On the other hand, the experimental set up described in the previous references 
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could be used for the realization of the QKR proposed in this work. The most 
important difference with the realized experiments would reside in that the 
intensity of the laser that hits the cloud of atoms should depend on time. 

In summary, we have developed a variant of the QKR model that shows 
anomalous diffusion in resonances. The spreading of the wave function has 
a rich variety of behaviors that go from the localized to the super-diffusive. 
This feature opens interesting possibilities for quantum information process- 
ing because it could be used to control the spreading of the wave function. For 
the primary resonances, in the long-time asymptotic regime, we have found 
an analytical expression for the dependence of the standard deviation with 
the strength parameter K{t). This expression shows explicitly how to select 
the unitary evolution to obtain a predetermined asymptotic behavior for the 
wave function spreading: super-ballistic, ballistic, sub-ballistic, diffusive, sub- 
diffusive, logarithmic or Bessel-localized. The parameter a determines the de- 
gree of diffusivity of the system in the long-time limit. We have shown nu- 
merically that the standard deviation for the secondary resonances has the 
same qualitative behavior as that exhibited by the primary resonances. The 
probability distributions are similar in the case of localization for both types 
of resonances, but in the other cases the profiles are very different. 

Finally, it is important to mention that there exists a parallelism between the 
results obtained in the present work and those obtained in Ref. [T7] for the 
QW. This evidences the theoretical equivalence between the QKR in resonance 
and the QW [71I2T1I20] . 
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